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Abstract 

We study the density of quasiparticle states of dirty ci-wave superconduc- 
tors. We show the existence of singular corrections to the density of states 
due to quantum interference effects. We then argue that the density of states 
actually vanishes in the localized phase as \E\ or E 2 depending on whether 
time reversal is a good symmetry or not. We verify this result for systems 
without time reversal symmetry in one dimension using supersymmetry tech- 
niques. This simple, instructive calculation also provides the exact universal 
scaling function for the density of states for the crossover from ballistic to 
localized behaviour in one dimension. Above two dimensions, we argue that 
in contrast to the conventional Anderson localization transition, the density 
of states has critical singularities which we calculate in a 2 + e expansion. We 
discuss consequences of our results for various experiments on dirty high-T c 
materials. 

I. INTRODUCTION 

The question of the quasiparticle density of states of a two dimensional d x 2_ y 2 supercon- 
ductor in the presence of disorder has been a matter of some controversy. Early theoretical 
work based on approximate self-consistent treatments of the disorder demonstrate that 



a finite Fermi level density of states is generated for arbitrarily weak disorder. In contrast, 
some exact results for a simplified model of the disorder which ignores the scattering 
between the two pairs of antipodal nodal points show that the density of states van- 
ishes on approaching zero energy (measured from the Fermi energy) as p(E) ~ E? . Claims 
of rigourous proofs of a constant non-zero density of states have also appeared in the 
literature. 

In a recent paper || , we discussed the problem of quasiparticle transport and localization 
in dirty superconductors ignoring the quasiparticle interactions, and treating the disorder 
with a non-linear sigma model field theory. The starting point for the sigma model descrip- 
tion is the approximate self-consistent treatment of the disorder which, as mentioned above, 
generates a finite density of states. We argued that inclusion of small harmonic fluctuations 
about the self-consistent solution leads to diffusion of the spin and energy densities of the 
quasiparticles (though not of the charge density). Quantum interference effects finally lead 
to quasiparticle localization at the longest length scales in two dimensions. In this paper, 
we consider the behaviour of the density of states in the sigma model. We show that in 
the diffusive regime, quantum interference effects lead to a singular logarithmic suppression 
of the density of states. We then argue that in the localized spin insulator, the density of 
states actually vanishes as \E\ for superconductors with both spin rotation and time rever- 
sal symmetry. A schematic plot of the density of states as a function of energy is shown in 
Figure |I[ The linear density of states of the pure d x 2_ y 2 superconductor gets rounded off at 
an energy scale E\ of the order of the elastic scattering rate. This marks the crossover from 
the ballistic to the diffusive regime. At a lower energy scale E 2 ~ p, the density of states 
dips linearly to zero. (Here D is the "bare" spin diffusion constant, and £ is the quasiparticle 
localization length). This second energy scale marks the crossover from the diffusive to the 
localized regime. The ratio of the two crossover scales jj± is exponentially large in the bare 
dimensionless spin conductance, and can be quite large. 
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FIG. 1. Density of states of the two dimensional dirty d x i_ y i superconductor 

Note that our conclusion of the vanishing density of states is only superficially similar 
to the results of Nersesyan et. al. [|J. In particular, we argue that the localization length 
approaches a finite constant as E — > 0. (In contrast, Nersesyan et. al. find a diverging 
localization length as E —>■ 0). 

For a superconductor with spin rotation invariance but no time reversal (T), there is 
again a logarithmic suppression of the density of states in the diffusive regime. This accounts 
entirely for the leading logarithmic suppression of the spin conductance found earlier, and 
provides an explanation of it. In this case, we argue that in the localized phase, the density 
of states vanishes as E 2 . 

We provide an explicit verification of some of our general results by exact non- 
perturbative calculations in one dimension using supersymmetry techniques. 
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II. DIFFUSIVE REGIME 



Consider first a two dimensional d x i_ y i superconductor with both spin rotation and T 
invariance. In the absence of quasiparticle interactions, and at scales larger than the elastic 
mean free path, the quasiparticle diffusion M is described by the replica field theory M 

S^LaM = J d 2 x^-Tr (W • Vf/ f ) - V Tr (U + rf) (1) 

Here U(x) are In x In unitary matrices G Sp(2n), and n is a positive infinetisimal that is 
introduced to enable calculation of the appropriate Green's functions. The coupling constant 
g is related to the spin conductance a s by - = |cr s . 

The quasiparticle density of states at the Fermi energy is exactly proportional to the 
uniform spin susceptibility which is the "order parameter" || for this field theory. Thus it 
is given by 

where po is the bare d.o.s (i.e, it's value on the scale of the mean free path). The limit 
n —>■ is implied. As this is the "order parameter" for the field theory, quantum corrections 
to it can be obtained from the known results || for the "field renormalization" to one loop 
order. To perturbatively calculate corrections to p, we write U = 1 + i(j) — %r + ■ ■ • with 
belonging to the Lie algebra of Sp(2n), and expand in powers of 0. To quadratic order, the 
action is 

S ° = kg I d2xTr ^^ ( 3 ) 

We may choose a basis T a for the Lie algebra to write <fi = 4> a T a , and normalize the basis 
by some choice Tr{T a T h ) = C5 ab with C a positive constant. The action then becomes 

S * = Y g j ■ V ^ ( 4 ) 

(Summation over the index a is implied in the above equations). The matrices T a are 
traceless, and the susceptibility may be expressed to leading order as 
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= ^{An-Tr(T a T b )<j )a( j )b } 

= 1-^G(0) 

In the last line, we have used the function G(x) defined by (0 a (x)0b(O)) = G(x)5 a b] N is the 
number of linearly independent matrices T a . For Sp(2n), N = n(2n + 1). It is understood 
that the n = limit is taken at the end of the calculation. Considering now a finite system 
size, we get 

r(o) = I [ d2k 1 

This gives 

P ~ P0 - 9N ln(±) (5) 
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Note that the constant C has dropped out of this result (as it should). The quantity l e is 
the elastic mean free path. 

The same result is obtained in the T broken but spin rotation invariant case. This is 
described by the s ffi) fi^d theory with the action: 

S = J d 2 x^Tr ((VQ) 2 - V Qa z ) (6) 

where Q = U^a z U with U G Sp(2n). The density of states is again the order parameter of 
this field theory, and is given by 

Calculation similar to the one above gives the result Eqn. [| but with N = n(2n + 1) — n 2 
being the number of independent massless fields. In either case, in the replica limit, we get 

P - P " 9 \4t) (8) 
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Thus to leading order, the suppression of the density of states is independent of whether or 
not T is present. 



The leading logarithmic correction to the spin conductance in two dimensions was eval- 
uated in Ref. 101: 



<7„ 



— n In T~ 

2tt 2 L 



(9) 



where cr° is the bare spin conductance, and l e is the elastic mean free path. If T is broken, 
then the correction is reduced by a factor of two H. 

The spin conductance satisfies the Einstein relation cr s = ^ with £) being the spin 
diffusion constant. (The factor of four arises from the spin of 1/2.) Eqns. (^) and (£5|) 
together with the relation ~ = |cr s imply a logarithmic suppression of the diffusion constant 
at order ~ when T is present. Without T invariance, there is no suppression of the diffusion 
constant to this order. 

It is possible to understand this result in terms of a semiclassical picture f£,10] involving 
interfering trajectories. To that end, consider, quite generally, a lattice Hamiltonian for the 
quasiparticles in a singlet superconductor 

^2 C l<T C j° + ^ij c l] c ]i + ^ij C ii C *T 
cr 

where i, j refer to the sites of some lattice. Hermiticity implies ty = and spin rotation 
invariance requires Ay = A^. It is useful conceptually to use the alternate representation 
in terms of a new set of <i-operators defined by: = c^, dii = c\<. The Hamiltonian Eqn.[T0| 
then takes the form 

f > \ 



(10) 



u 



tin ^-ii 



Writing „ tJ 
antisymmetric, we get 



ij 

ia v i4 with real symmetric and bij 



(11) 



A* +* 
\ ^ij ''J ) 



Hij = iby + dlj.a 



-bji, real 



(12) 



Note that SU(2) invariance requires a y Hija y = —H*y This implies that the amplitude 
iGij,ap = (ia\e~ lHt \j (3)9(t) for a d-particle to go from point j, (pseudo)spin (3 to point i, 
spin a satisfies the relations: 
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GW*) = -G* m (t) (13) 
<W*) =G*^(t) (14) 



The Fourier transform of this amplitude is 



G ijtaP (u + iri) = / dte^+^G^t) 



= FT i • lift) 

w — -H + zr/ 

The density of states at the Fermi energy may be obtained from this in the usual manner. 

p = _I/ m (^ n (^) + (T~j) (15) 

Consider now the return amplitude Gu^(t). This can be written as a sum over all 
possible paths for this event. Consider in particular the contribution from the special class 
of paths where the particle traverses some orbit and returns to the point % in time t/2 with 
spin down, and then traverses the same orbit again in the remaining time and returns with 
spin up. This contribution to iGa^{t) can be written 



iG ii,U (2) iG Wt - - \ G ii,ll (2 



|2 



using the symmetry relation Eqn.|T4[ Now iG^jd)! 2 is just the probability for the event 
i f— > % I in time t/2. For large t, this is half the total return probability which ~ | in 
two dimensions if the particles are diffusing. This leads to a logarithmic divergence in the 
density of states which may be cutoff by a finite system size. To be precise, this gives 

yr-^ ln (i) (16) 

in agreement with the field theoretic result obtained earlier. 

In addition, even to leading order the spin conductance for T invariant systems is sup- 
pressed further by the usual constructive interference between paths and their time reverse 
which explains the larger suppression in that case. 



7 



III. LOCALIZED REGIME 



Having established the presence of a singular suppression of the density of states in 
perturbation theory, we now consider the opposite limit of strong disorder when the system is 
localized. We show that the density of states vanishes at zero energy. To see this heuristically, 
consider the Hamiltonian (PD| ) in the limit of strong on-site randomness and weak hopping 
between sites. In the limit of zero hopping, the sites are all decoupled. At each site, the 
Hamiltonian in terms of the ci-particles satisfies the SU(2) invariance requirement a y Ha y = 
—H*. This takes the form H = a. a with a random. With T symmetry, we further have 
H = H* implying a y = 0. Considering now the case where the probability distribution of a 
has finite, non-zero weight at zero, we see immediately that the disorder averaged density of 
states vanishes as E 2 without T and as \E\ with T. Now consider weak non-zero hopping. 
In the localized phase, perturbation theory in the hopping strength should converge, and we 
expect to recover the single site results at asymptotically low energies. 

A more formal field theoretic version of this argument with the same conclusions is as 
follows. As we are concerned with the properties of the localized phase, we prefer to phrase 
the argument in terms of a supersymmetric field theory rather than the replica version used 
before. In the localized phase, we expect that a strong coupling expansion of this field 
theory converges. This may be performed, as usual, by regularizing the sigma model on 
a lattice. The leading term in the strong coupling expansion is the zero dimensional limit 
of the sigma model which is equivalent to the random matrix theory of Hamiltonians with 
these symmetries. In the random matrix limit, it is known that the density of states 
vanishes in the manner disussed above. 

These results on the localized phase can be verified in great detail in one spatial dimension 
for systems without T. Consider a lattice Hamiltonian for the ci-particles in one dimension. 
In the absence of disorder, we take this to be of the form 

7i = — t(d\a z d i+ i + h.c) — fid\a z di (17) 
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Now consider adding random terms to this Hamiltonian consistent with the required 
symmetry. For weak disorder, we may just keep the modes near the two Fermi points of the 
pure system. Linearizing the dispersion near these Fermi points, we arrive, as usual, at a 
one dimensional Dirac theory with various sorts of randomness. The resulting Hamiltonian 
can be written down on symmetry grounds as 

H = -ir z d x + (fh(x).ff) (Ap-) ~ {V2{x).a) + t (x)T y + (t(x).a)r x (18) 

This is the most general Hamiltonian consistent with the symmetry a y Ha y = — H* required 
by spin rotation invariance. (We have set the Fermi velocity to one). The f are Pauli matrices 
in the right mover/left mover space and r]i,r] 2 ,to,t are random, independently distributed 
real variables. Green's functions of this Hamiltonian are generated by the action 



S = J dx [ip(iH + + C{iH + (19) 



where ip, ip are Grassmann variables, and £ is a complex scalar field. For a system of 
finite size L, we impose periodic boundary conditions on all fields. The partition function 
Z corresponding to this action is exactly equal to one for any L as the fermionic and 
bosonic integrals cancel each other. Here u is chosen to have a positive real part to ensure 
convergence of the bosonic integral. 

The density of states can be obtained from the Green's function through 



p(E) = Im TrGiE + irj) (20) 

7T 

where the overline indicates disorder averaging and 

G${x^E + irj) = (aax\ 1 \bf3x') (21) 

hj — ri + 17] 

It's disorder average can be expressed in terms of correlators of the either the Bose or Fermi 
variables: 



iG a J(x,x';E + i V ) = (^ aa {x)^Ux')) = (U(x)C b p(x')} (22) 



9 



(We have set iuo = E + irj in evaluating the correlators). As we need the Green's function 
when x = x', there is some subtlety on the relative ordering of x and x'. The correct 



procedure [jnj is to take a symmetrized form: 

2np(E) = Re [<£ OQ (x + e)£ a (x) + U(x ~ e)CW)l (23) 

where e = + , and summation over a, a is implied. Precisely the same expression with 
£ — > ip holds in terms of the fermionic variables as well. The other physical quantity we 
will be interested in is the diffusion propagator. This is defined , as usual, in terms of the 
Green's function by 



P(x,x') = £ \G aa;b p(x,x';iri)\2 (24) 

ab,afi 

Now the symmetry a y Ha y = —H* can be used to show that 

G:f(x,x';i V ) = -(-ir^G:i(x,x';i V ) (25) 
where a = 2 if a = 1 and vice versa. Thus P(x, x') may be written 



P(x,x') = -j:(-ir^G a jG^ (26) 

ab,a(3 

= E (-l) a+/3 (^ a (x)^(a;')e ac5 (a;)^(a;')) (27) 

ab,a(3 

We have chosen to write one Green's function in terms of the fermions and one in terms 
of the bosons. This enables a calculation of the two particle properties using the same 
formulation needed to calculate the one particle properties. 

In the limit where f 1 = 0,/i = 0, 1,2,3, the left/right moving fields decouple for every 
realization of the disorder. Considering just one of them, say the right-movers, we get the 
action 



S = dx 



ij 1 {d x + fj 1 .a)^ l + u{^ l ) + ^j ^i) (28) 



We now average over the disorder assuming 771 to be distributed as P[fji] oc exp — J dx 



Am) 2 



2u 



The resulting translationally invariant action can be interpreted as the coherent state 
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path integral of a zero dimensional quantum "Hamiltonian" in terms of bose operators 
bia = (&it>&ij.) an d fermi operators f\ a = Before doing that, we note that 

the fermionic fields actually satisfy periodic boundary conditions. To get fermion fields 
that satisfy antiperiodic boundary conditions, we may perform a change of variables 
4>i - ► ^Pi^ 1 ^ , fpi V'le -1 ^. This adds a term / dx^ipitp to the action. Thus, we get 

Z = STre- LhR (29) 
h R = u (fla^ + (/i <- b,)) 2 + utflh + 6}6i) (30) 

(The subscript R on /i is a reminder that this is for the right-moving fields alone). The 
supertrace operation STr is defined through STrO = Tr ((— l)-^ 1 ^. It is necessary to 
take the supertrace to account for the extra term in the action coming from the change of 
the fermion boundary conditions. 

At zero u, it is clear that there is a triplet of states with zero energy: the vacuum 
state with no particles which we denote |0), the state /if/ijO) = | Tl;0), and the state 
^2(^11 fli — biifii)\0) = I |; I)- All other states have energies at least 0(u). (Non-zero uo 
of course splits the energies of this zero energy triplet). Similar considerations apply to the 
left moving sector as well. Thus there are a set of nine states all at zero energy at zero 00 in 
the limit of decoupled right/left sectors. 

Now consider coupling the left/right moving sectors. The full action can also be inter- 
preted (after disorder averaging) as the coherent state path integral of a zero dimensional 
quantum Hamiltonian. In the limit where the coupling is small, it is sufficient to project the 
interactions induced between the two sectors to the nine dimensional space of the ground 
states of the two decoupled sectors. For simplicity, we assume that the are Gaussian 
distributed with 

t^x)t u (x') = t 2 5^5(x - x') (31) 

To leading order then, the coupling between the two sites in the nine dimensional space will 
be of order ^. To derive the form of this coupling, it is convenient to gauge away ffi and 
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for 



772 by letting ^ -> L/^,^ -> V^;, and similarly for £ with ^(x) = T x [ e */«k*0<0-* 
i — 1,2. (Tj. is the x-ordering symbol.) We impose the condition that Ui(x — L) — 1 to 
maintain the periodic boundary conditions. Note that the Ui(x) are random 5*f/ (2) matrices. 
The full action can then be written 



S= I dx ij){T z d x + + i{^B{x)T + iP + ij)B\x)T-ij)) + -»• (32) 
Here B(x) = Ui{x){t + it- a)U\ (x) is a random 2x2 matrix. It is distributed according to 



B a0 (x)B l5 (x') = t 2 e- ulx - x ' l {a y ) aS {a y )^ 

For large u, we may replace t 2 e - u \ x - x '\ _ > J5(x — x') with J — It is now convenient to 
change variables ip2a ~i J 2ai i , 2a ~^ tyim £20 — Go- This changes the action to 

| [#> 2 v + r^,e + (-TiB(xW2 - ^(x)^ - + CsB^i) + w(# + ro] (33) 

Under this change of variables, the expression Eqn. |23] for the density of states remains 
unchanged (in the limit e — > + ). 

We may now perform the disorder average to get a translationally invariant action which 
can be interpreted as the coherent state path-integral of a zero dimensional quantum problem 
with bose operators b aa and fermi operators / aQ (a = 1, 2; a =j, j) and a "Hamiltonian" : 

h = h + h u) (34) 

to = -J {[{fl°yfl){ft°yft) + (/l*y/l)(M,/ 2 ) + 2(A t / 1 - 1)(/ 2 V 2 - 1); 
-2(61^ + 1)(6 2 6 2 + 1) + [2{b\a y fl)Uia y bl) + 2(b 1 a y f 1 )(f 2 a y b 2 ) 
-2(/ 1 t 6 1 )(/ 2 t 6 2 )-2(6 t J 1 )(6 2 / 2 ))]} (35) 

^ = ^(/ t / + & t &) (36) 

Note that this Hamiltonian is non-hermitian. It's action on the nine dimensional subspace 
is simplified by noting that f\fi — f 2 f 2 , an d b\bi — b 2 b 2 commute with h. Thus the six states 
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I T;i>i®|0) 2 ,|0)i<8>| T;l>2,| U;0)i®|o) 2 ,|o)i(8)| U;0) 2 ,| T;l)i®| U;0) 2 ,| U;0)i®| T;|) 2 

are immediately seen to be eigenstates of the Hamiltonian. The first four have eigenvalues 
4 J + 2uj and the last two have 4 J + 4u. The action of h on the remaining three states 
1 0) i <g> |0) 2 , | T;!)i ® | T; D2, | U;0)i <g> | Tl;0) 2 can be represented in terms of a 3 x 3 



non-hermitian matrix: 



\ 



V 



(37) 



/ 



-J -J 
J 2J + uj J 
—J —J UJ 

There is one eigenvalue and two eigenvalues 4 J + Au. With these eigenvalues, it is easy 
to see that Z = STre~ Lh = 1 for any system size L, as required. We will also need to know 
the zero energy wavefunction. However due to the non-hermiticity of the Hamiltonian, the 
left eigenvector (Vl) an d the right eigenvector (Vr) are different. They are easily seen to be 



V R 
Vl 



1 + 



1 



i + £ Li 



We have normalized these so that Vl-Vr = 1. There are some subtle questions regarding the 
resolution of the identity in the basis of right (left) eigenvectors of the Hamiltonian which 
are addressed at length in the Appendix. 

The expression Eqn. for the density of states can clearly be interpreted as the following 
expectation value: 



2vrp = 2Re(2 + 6 f 6) 



(38) 



where we calculate expectation values setting i& = E+irj. In the thermodynamic limit only 
the zero energy state contributes, and the result is 



1 _ & 

2^ } (1 + f) 2 



(39) 



Note that this vanishes as E 2 at small E, entirely consistent with the general arguments 
presented earlier. For large E, this saturates at - which is the ballistic result (see Figure 0). 
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FIG. 2. Density of states in one dimension in the absence of time reversal symmetry; the energy 
scale J ~ | 

The diffusion propagator can also be calculated explicitly in this one dimensional case. 
The calculation proceeds straightforwardly from Eqn. |26|. We first perform the change of 
variables ip 2a — > ~i J 2a i ~~ > ^la^ia — CSa; anc ^ then interpret the resulting correlator 
as an expectation value of an operator in the equivalent quantum problem. We find 

P(x,x') = (T x (0 1 (x)0 2 (x'))) 

Oi = (fnhi - fiM + (f 2 / 2l - fiA) 

The expectation value is to be taken in the zero energy state. Consider x > x' for definiteness. 
Thus, we write 

P(x,x') = (0 1 e- h(x ~ x ' ) 2 ) (40) 

We may evaluate this by inserting a complete set of states. As 2 acting on the ground 
state is a state with energy AJ (when u — > 0), P(x,x') decays as e~ 4J< ^ x ~ x '\ The precise 
result is easily seen to be 
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P(x, x) = 8e~ iJlx - x ' 1 (41) 
Thus the localization length of the system is £ = jy. In momentum space, this becomes 

™ = (4^ ^ 

Note the difference in structure from conventional localization with a finite density of states 
where P(q = 0,u) has a pole at u = 0. In this problem, the density of states vanishes and 
there is no pole. 

The detailed calculation above of the one dimensional problem is strong evidence in 
support of our general assertions regarding the vanishing of the density of states in the 
localized phase. In this case, the crossover from the constant to the vanishing density of 
states occurs at an energy scale J ~ | which is the energy scale for the crossover from the 
ballistic to the localized regime. We expect that Eqn. [39] is a universal scaling function 
for the density of states associated with this crossover. In two dimensions (or in quasi one 
dimensional situations such as that considered in Ref. flI2|), the crossover occurs between 
the diffusive and localized regimes at a scale p (where D is the diffusion constant). Again, 
this crossover is expected to be represented by a universal scaling function for the density 
of states. 



IV. ABOVE TWO DIMENSIONS 

We now turn to the situation above two dimensions where there is the possibility of a spin 
metal to spin insulator transition. The density of states is finite in the spin metal phase and 
vanishes on approaching the transition. Thus, in contrast to usual Anderson localization, 
the density of states behaves as a conventional order parameter in these universality classes. 
The order parameter exponent (3 may be calculated within the 2 + e expansion. We find, 
to leading order in e, (3 — ~ if T is present, and (3 = 1 without T. Right at the transition, 
the density of states vanishes with energy as p{E) ~ E~$ . The exponent 8 = - with and 
without T respectively. 
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V. DISCUSSION 



In this paper, we have studied the behaviour of the quasiparticle density of states in a 
dirty d x 2_ y 2 superconductor ignoring the quasiparticle interactions. We showed the existence 
of a singular logarithmic suppression of the density of states in the diffusive regime in two 
dimensions due to quantum interference effects. We then argued that in any dimension 
in the localized phase the density of states vanishes as \E\ if both spin rotation and T 
symmetry are present, and as E 2 if spin rotation is the only symmetry. This was verified 
by a simple explicit calculation in the latter case in one dimension using supersymmetry 
techniques. Above two dimensions, we showed that the density of states is finite in the spin 
metal phase, but vanishes on approaching the transition to the insulator. The corresponding 
critical exponent was calculated in a 2 + e expansion. These results are summarized in Table 

I 

These results imply that the spin susceptibility, linear temperature coefficient of the spe- 
cific heat, and the tunneling density of states all have a logarithmic suppression as a function 
of temperature in the diffusive regime in two dimensions due to quantum interference effects. 
As pointed out in Ref. ||, inclusion of a Zeeman magnetic coupling drives the system into 
the usual unitary universality class where there are no singular corrections to the density of 
states. Thus this logarithmic correction is killed by an external Zeeman field (though not 
by a purely orbital magnetic field). Experimental verification of this effect may be clouded 



somewhat due to the presence of quasiparticle interactions. We have shown elsewhere [13 
that in the diffusive regime, interaction effects lead to a logarithmic Altshuler-Aronov sup- 
pression of the tunneling density of states in the diffusive regime in two dimensions. This 
therefore adds to the quantum interference correction discussed in this paper. In contrast, 
the specific heat and spin susceptibility are expected to get logarithmic enhancements due 



to interactions in two dimensions in the diffusive regime ||14| , which too is killed by a Zeeman 
field. They would thus compete with the quantum interference corrections. Nevertheless, if 
the interactions are weak, we expect that the quantum interference effects would dominate 
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leading to a logarithmic suppression of the spin susceptibility and specific heat, which can 
be probed by applying an external Zeeman magnetic field. 

The effect of interactions in the localized phase is a more delicate matter. Qualitatively, 
repulsive interactions tend to favor the formation of local moments leading possibly to a 
divergent spin susceptibility and linear specific heat coefficient. This effect will however 
compete with the vanishing density of states we have discussed above (which tends to pro- 
duce a vanishing spin susceptibility, etc). The ultimate fate of the localized phase in the 
presence of these two competing physical effects is a formidable problem that we will not 
attempt to answer here. 

We are particularly grateful to Martin Zirnbauer for a most useful communication. We 
also thank Leon Balents, Ilya Gruzberg and Andreas Ludwig for useful discussions. This 
research was supported by NSF Grants DMR-97-04005, DMR95-28578 and PHY94-07194. 
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TABLES 



TABLE I. Properties of the two different symmetry classes of superconductors with spin SU(2) 
symmetry. WL stands for weak localization. pi oc {E) is the density of states in the localized phase. 
The last column gives the critical properties of the density of states above two dimensions as 
calculated in a 2 + e expansion. The distance from the critical point is (Sg), and p cr (E) is the 
density of states at the critical point. 



SYMMETRY 


WL in d = 2 


Ploc(E) 


CRITICAL PROPERTIES IN d = 2 + e 


Spin SU{2) and T 


Sp _ 1 i„ ( L\ 
p 4n' 2 a s m \l e ) 




\ E \ 


p(£ = 0)~ ($0)5 

p cr (£) ~ i^i f 


Spin SU{2) and no T 


§£- i l n f^ 

p 4tt 2 o- s \l e J 




E 2 


p(£ = 0) ~ (<5 5 ) 

Pcr(E) ~ |£|§ 
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APPENDIX A: RESOLUTION OF THE IDENTITY 



In this appendix, we discuss some subtle questions regarding the resolution of the identity 
in the eigenbasis of the (super) Hamiltonian Eqn. [33]. As all the subtleties are associated with 
the three dimensional subspace spanned by |0)i <S> |0) 2 , | f ; |)i ® | |; j) 2 , | |J.; 0)i <S> | ||; 0) 2 , 
we just focus on these three states. In this subspace, the Hamiltonian h is represented by 
the 3x3 matrix Eqn. |37|. The right eigenstates corresponding to the two eigenvalues and 
4 J + 4u; are easily seen to be (in bra/ket notation) 



|-Ri) = ai 



l + z 
-1 
1 







(Al) 



where ai,a 2 are normalization constants, and z = j. {\Ri) has eigenvalue 0, and |i? 2 ) has 
eigenvalue 4 J + Au). Note that the right eigenstates do not span the full three dimensional 
space. From the structure of the Hamiltonian Eqn.|34], it is easy to see that from every right 
eigenstate \R), a corresponding left eigenstate (L\ can be obtained by the operation 



(L| = ((-1)4*|J2» 5 



(A2) 



where the symbol T denotes taking the transpose. The left eigenstates corresponding to 
\Ri) and \R 2 ) then are 



(L 1 \=a 1 [l + z, 1, 1]; (L 2 \=a 2 [0, 1, 1] 



(A3) 



respectively, as can also be seen by direct calculation. The left eigenstates also do not span 
the full three dimensional space. Note that (Li\R 2 ) = (L 2 \Ri) = (L 2 \R 2 ) = 0. 

To get a complete set of states, we need to supplement (Li\ and (L 2 \ by any other linearly 
independent bra vector (La|. It is convenient to choose this to be orthogonal to \R\) and 
{{L 2 \f. 



(L 3 \ = a 3 [-2, -(l + z), (l + z)] 



(A4) 
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A corresponding right state \Rg) can be defined using Eqn. 



I-R3) = «3 



-2 
l + z 
l + z 



(AS) 



Clearly the \Ri){i = 1,2,3) form a complete set of states (as do the (Li\). By construction, 
we have the relations (L 3 |i?i) = (Li\R 3 ) = 0. We now impose the normalization conditions 



(LtlRx) = (L 2 \R 3 ) = (L 3 \R 2 ) = (L 3 \R 3 ) = 1. This fixes 01 
now possible to construct the resolution of the identity 



1+2' 



a 2 



l+z' 



0,3 



1 = |i*i)(Li| + |i22)(L 3 | + 1^3)^21 - |i? 2 )(i 2 



i It is 



(A6) 



This can be checked directly by it's action on any vector in the three dimensional space. 

This resolution of the identity can now be used to easily show explicitly that Z = 
Stre~ hL = 1 in the full nine dimensional space. For the calculation of the density of states 
or the diffuson, we need to know the action of the Hamiltonian h on I-R3). This is easily 
seen to be 



h\R 3 ) = U(l + z)(\R 3 ) + \R 2 )) 



Combined with the eigenvalue equation h\R 2 ) = 4/(1 + z)\R 2 ), this implies that 



(A7) 



e- hL \R 3 ) = e- AJ ^ L {\R 3 )-AJLz{l + z)\R 2 )) 

In the limit L — > 00, e~ hL \R 3 ) — > 0; similar considerations apply to (L 3 \ as well. Calculation 
of any correlation function is thus reduced, in the limit of infinite system size to a calculation 
in the zero energy state with (right) eigenvector \Ri). 
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